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ABSTRACT 

We explore the cosmological consequences of interacting dark energy (IDE) models using the SNLS3 supernova samples. In partic¬ 
ular, we focus on the impacts of different SNLS3 light-curve fitters (LCF) (corresponding to “SALT2”, “SiFTO”, and “Combined” 
sample). Firstly, making use of the three SNLS3 data sets, as well as the Planck distance priors data and the galaxy clustering data, 
we constrain the parameter spaces of three IDE models. Then, we study the cosmic evolutions of Hubble parameter H(z), decelera¬ 
tion diagram q(z), statefinder hierarchy S^\z) and S^\z), and check whether or not these dark energy diagnosis can distinguish the 
differences among the results of different SNLS3 LCF. At last, we perform high redshift cosmic age test using three old high redshift 
objects (OHRO), and explore the fate of the Universe. We find that, the impacts of different SNLS3 LCF are rather small, and can not 
be distinguished by using H(z), q(z), S^\z), S^\z), and the age data of OHRO. In addition, we infer, from the current observations, 
how far we are from a cosmic doomsday in the worst case, and find that the “Combined” sample always gives the largest 2cr lower 
limit of the time interval between “big rip” and today, while the results given by the “SALT2” and the “SiFTO” sample are close to 
each other. These conclusions are insensitive to a specific form of dark sector interaction. Our method can be used to distinguish the 
differences among various cosmological observations. 

Key words, cosmology: dark energy, observations, cosmological parameters 
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1. Introduction 

Type la supernova is a sub-category of cataclysmic variable stars 
that results from the violent explosion of a white dwarf star in a 
binary system ( Hillebrandt, & Niemeyerf2000[ ). Now it has be¬ 
come one of the most powerful probes to illuminate the mystery 
of cosmic acceleration ( |Riess et al.|1998[|Perlmutter et al.| 1999| >, 
which may be due to an unknown energy component that can 
produce an anti-gravitational effect, i.e., dark energy (DE), or a 
modification of general relativity, i.e., modified gravity (MG). 
[[] Along with the rapid progress of supernova (SN) cosmology, 
several high quality SN datasets were released in recent years, 
such as “SNLS” jAstier et~aL ( 2006|l, “Union” [Kowalski et al. 


et al. 


Suzu ki et al.| ( [20T2] >, and “Pan-STARRSl” |Scolnic et al.| (|2014). 


( 20081 “Co nstitution” |Hicken et al.| (|2009l), “SP SS” |Kessler 


( 2009]) “Uni on2’’[Amanull ah et al.|(|2010), “Union2.1 ” 


In 2010, the Supernova Legacy Survey (SNLS) group released 
their three years data |Guy et ah] ( |2010[ ). Soon after, combining 
these SN samples with various low-z to mid-z samples and mak¬ 
ing use of three different light-curve fitters (LCF)pj Conley et al. 
presented three SNLS3 data sets |Conley et al.[f2011) >: “SALT2”, 
which consists of 473 SN; “SiFTO”, which consists of 468 SN; 


1 For recent reviews, see Frieman, Turner, & Huterer 

2008);|Caidwell 

& Kamionkowski (2009); Uzan (2010); 

Wang (2010); 

Li, Li, Wang, & 

Wang (2011 2013); Weinberg et al. (20 

*-2 — % r 'nn “ Jh-rrU-r-rr a* At 

13). 


and “Combined”, which consists of 472 SN. It should be men¬ 
tioned that, in the cosmology fits, the SNLS group treated two 
important quantities, stretch-luminosity parameter a and color- 
luminosity parameter /3 of SNe la, as free model parameters. No¬ 
tice that a and are parameters for the luminosity standardiza¬ 
tion. In addition, the intrinsic scatter <x,„ f are fixed to ensure that 
X 1 Ido f = 1. 

A most critical challenge of SN cosmology is the control of 
the systematic uncertainties of type la supernovae (SNe la). One 
of the most important factors is the potential SN evolution. For 
examples, previous studies on the Union2.1 Mohlabeng & Ral- 
|ston| ( |2013] ) and the Pan-STARRS 1 data sets |Scolnic et al.| ( ]2014[ i 
all indicated that /? should evolve along with redshift z. Besides, 
it was found that the intrinsic scatter cr int has the hint of redshift- 
dependence that will significantly affect the results of cosmology 
fits ( |Marriner et al.|20TT| l. In addition to the SN evolution, an¬ 
other important factor is the choice of LCF (Kessler et al.|2009]). 


For instance, it has been proved that even for the same SN sam- 
plesjrsing “MLCS2k2” ( )Jha et al.fcOOTj l and “SALT2” ( jGuy etj 


various cosmological models (Sollerman et al. 

2009; Sanchez 

et al. 2009 

Pigozzo et al. 2011, Smale & Wilts 

lire 2011 Ben- 

gocheapO 

JBengochea & De Rossi|(2014 ). 


can estimate the distance information of each SN. 


One of the present authors has also done a series of research 
works to study the systematic uncertainties of SNe la. By using 
the SNLS3 dataset, we found that a is still consistent with a con¬ 
stant, but /? evolves along with z at very high confidence level 
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(CL) |Wang & Wang| ( (2013a)) . Soon after, we showed that this 
result has significant effects on the parameter estimation of stan¬ 
dard cosmology ( Wang, Li & Zhang|2014| >, and the introduction 
of a time-varying /3 can reduce the tension between SNe la and 
other cosmological observations. Moreover, we proved that our 
conclusion holds true for various DE and MG models (|Wang, et| 


|al.(2014j [Wang, Wang & Zhang|2014l|Wang, et aLfcOBj ). Be 

sides, the evolution of a has also been found for the JLA sampls 
(Li et al. 20161. In addition, by using three different SNLS3 
LCF, including “SALT2”, “SiFTO” and “Combined”, we briefly 
discussed the effects of different LCF on the parameter estima¬ 
tion of the A-cold-dark-matter (ACDM) model Wang & Wang 


(201 3a(> and the holographic dark energy (HDE) model (|Wang, 


et al. |2015j >. It must be emphasized that in these two papers 


only the fitting results given by different SNLS3 LCF are shown, 
and the corresponding cosmological consequences are not dis¬ 
cussed. Therefore, the impacts of different LCF are not studied 
in details in our previous works. The main scientific objective 
of the current work is presenting a comprehensive and system¬ 
atic investigation on the impacts of different SNLS3 LCF. To do 
this, both the cosmology-fit results and the cosmological conse¬ 
quence results are shown in this work. 

As mentioned above, the effects of different SNLS3 LCF on 
the ACDM model has been briefly discussed in (Wang & Wang 
(2013a]) . To obtain new scientific results, new elements need to 
be taken into account. Since the interaction between different 
components widely exist in nature, and the introduction of a in¬ 
teraction between DE and CDM can provide an intriguing mech¬ 
anism to solve the “cosmic coincidence problem” (|Guo, Ohta, & 


|Tsujikawa|2007||Li et al.|2009a|b He et al.|2010]) and alleviate 


the “cosmic age problem” (Wang & Zhang 2008; Wang, Li & Li 


2010||Cui & Zhang|2010 Duran & Pavon|201 1| >, it is interesting 

to study the effect of different LCF on interacting dark energy 
(IDE) model. Here we adopt the w-cold-dark-matter (wCDM) 
model with a direct non-gravitational interaction between dark 
sectors. It must be emphasized that, if our conclusions are de¬ 
pendent on the specific form of dark sector interaction, these 
conclusions will not be reliable at all. So in this work, three 
kinds of interaction terms are taken into account to ensure that 
our study is insensitive to a specific interaction form. In addi¬ 
tion, to make a comparison, we also consider the case of wCDM 
model without interaction term. 

According to the previous studies on the potential SN evo¬ 
lution, we adopt a constant a and a linear (5{z) — fio + [3\z in 
this work. Making use of the three SNLS3 data sets, as well as 
the latest Planck distance prior data ( |Wang & Wang|2013b| ), the 
galaxy clustering (GC) data extracted from Sloan digital sky sur¬ 
vey (SDSS) data release 7 (DR7) (Hemantha, Wang & Chuang 
|2014| ) and data release 9 (DR9) ( [Wang |2014j ), we constrain the 
parameter spaces of the wCDM model and the three IDE models, 
and investigate the impacts of different SNLS3 LCF on the cos¬ 
mology fits. Moreover, based on the fitting results, we study the 
possibility of distinguishing the impacts of different LCF by us¬ 
ing various DE diagnosis tools and cosmic age data, and discuss 
the possible fate of the Universe. 

We present our method in Section[2] our results in Section[3j 
and summarize and conclude in Section[4] 


2. Methodology 

In this section, firstly we review the theoretical framework of the 
IDE models, then we briefly describe the observational data used 
in the present work, and finally we introduce the background 
knowledge about DE diagnosis and cosmic age. 
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2.1. Theoretical Models 

This work is a sequel to the previous studies of our group, and 
thus we study the same IDE models considered in (Wang, et 
|ak||2014| ). It must be mentioned that, here we consider a flat 
Universe and treat the present fractional density of radiation as 
a model parameter, which is a little different from the case of 
(Wang, et al. 2014). In a flat Universal the Friedmann equation 
is 


3M~ pl H~ - Pc + Pde + Pr + Pb, (1) 

where H = d/a is the Hubble parameter, a = (1 + z) 1 is the 
scale factor of the Universe (we take today’s scale factor ao = 
1), the dot denotes the derivative with respect to cosmic time 
t, M p = (SnG) 1 is the reduced Planck mass, G is Newtonian 
gravitational constant, p c , p t j,., p r and pb are the energy densities 
of CDM, DE, radiation and baryon, respectively. The reduced 
Hubble parameter E(z) = H(z.)/Ilo satisfies 

E 2 = G c0 — + Ode o — + Oo — + fiw —, (2) 

PcO Pde 0 PrO PbO 

where Q c o, Tide o, Lf-o and Q/,o are the present fractional den¬ 
sities of CDM, DE, radiation and baryon, respectively. As 
mentioned above, we treat Q,o as a model parameter in this 
work. In addition, p r = p,o(l + z) 4 , Pb = pw>(l + z) 3 - Since 
Q.d e o = 1 - Lf d - G/xi - G,-o, Tide o is not an independent parame¬ 
ter. 

In an IDE scenario, the energy conservation equations of 
CDM and DE satisfy 

Pc + 3 Hp c = Q, (3) 

Pde + ?>H{p de + Pde) = ~Q, (4) 

where p c i e = wpd e is the pressure of DE, Q is the interaction term, 
which describes the energy transfer between CDM and DE. So 
far, the microscopic origin of interaction between dark sectors 
is still a puzzle. To study the issue of interaction, one needs to 
write down the possible forms of Q by hand. In this work we 
consider the following three cases: 


Q\ = 3 yHpc, 

Qi = lyHpde, 

1 u PcPde 

Qi = 3 yH- 


' Pde 


(5) 

(6) 
(7) 


where y is a dimensionless parameter describing the strength of 
interaction, y > 0 means that energy transfers from DE to CDM, 
y < 0 implies that energy transfers from CDM to DE, and y = 0 
denotes the case without dark sector interaction, i.e. the wCDM 
model. Notice that these three kinds of interaction form have 
been widely studied in the literature (Guo, Ohta, & Tsujikawa] 


20071 |Li et al.|2009b[|He et al.|2010||Ei& Zhang 201 4). Forsim- 

plicity, hereafter we call them IwCDMl model, I wCDM 2 model, 
and IwCDM3 model, respectively. 

For the wCDM model, we have 


E(z) = (Oro(l + z) 4 + (O w , + Q c0 )(l + z) 3 + fl*o( 1 + z) 3(1+M °)‘ /2 . 
For the IwCDMl model, the solutions of Eqs. ([3]) and (|4]i are 
Pc=p c od+z) 3(1 " r) , (9) 

3 The assumption of flatness is motivated by the inflation scenario. For 
a detailed discussion of the effects of spatial curvature, see jClarksoru] 
Cortes, & Bassett|2007| 
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Pde = -^((1 +Z) 3(1+W) - (1 +Z) 3(1 -^) + Pde0 ( 1 +Z) 3(1+H ' ) . (10) 
w + y 

Then we have 

E(z) = (Oo(l + z) 4 + £2w(l + z) 3 + £2*o(l + z) 3(1+H ° 

+Q c0 (—-—(1 + z) 3(1+w) + -^(1 +z) 3(1 ' r) )) 1/2 . (11) 

w + y w + y ' 

For the IwCDM2 model, the solutions of Eqs. Q and Q are 

P*=p rfe o(l+z) 3(1+H ’ +r) , (12) 


p c = p cQ ( 1 + z ) 3 + 1 + z f _ 2^(1 + Z )3(l+w+y) (13) 

w + y w + y 

Then we get 

E(z) = (Oo(l + z ) 4 + (Q f o + Q/,o)(l + z ) 3 

+n*o(^-(l + z ) 3 + -^-(1 + z) 3(1+w+ ^)) 1/2 . (14) 

w + y w + y ' 

For the IwCDM3 model, Eqs. 0 and ([4]) still have analytical 
solutions 

Pc — Pco(l + z) 3 ( Pc ° + Pde ° (1 +z) 3 (w+r) r^, (15) 

PeO + Pde 0 PcO + Pde 0 


Pde = Pde 0d +Z) 3(1+W+y) ( PC ° 

PcO + Pde 0 

+_ £^2. _(1 +z )3(H'+y))-^_ 

PcO + PdeO 

Then we obtain 

E(z) = (£2,0(1 + z) 4 + £1«)(1 + z) 3 + QcoC(z)(l + z) 3 
+Q* 0 C(z)(l + z) 3 a+w +r ^i/2_ 

where 


C(z) = ( 


a 


■cO 


a 


■de 0 


Q C 0 + Clde 0 £lcO + EldeO 


-(1 + z ) 3(m ' + ^)- 


(16) 


(17) 


(18) 


For each model, the expression of E(z) will be used to calculate 
the observational quantities appearing in the next subsection. 


2.2. Observational Data 


In this subsection, we describe the observational data. Here we 
use the SNe la, CMB, and GC data, which is same with our 
previous paper ( |Wang, Li & Zhang|2014) l. It should be empha¬ 
sized that, in this work we use all the three SNLS3 data sets (i.e. 
“SALT2”, “SiFTO”, and “Combined ” sam ple), while only the 
“combined” data are used in ( Wang, Li & Zhang[2014[ >. in addi¬ 
tion, the GC data are also updated in this work, compared to that 
used in ( |Wang, Li & Zhang|2014[ >. 


2.2.1. SNe la Data 


As mentioned above, we use all the three SNLS3 data sets. In 
the following, we briefly introduce how to include these three 
data sets into the^f 2 analysis. 

Adopting a constant a and a linear /J(z) = {3q + P\Z„ the pre¬ 
dicted magnitude of an SN becomes 

Wmod = 5 log 10 D l (z) - a(s - 1) +/3(z)C + M. (19) 

The luminosity distance Dpz) is defined as 

£>l(z) = H 0 ( 1 + ZhelMz), (20) 


where z and Zhei are the CMB restframe and heliocentric redshifts 
of SN, and 

■"-*£■& (21) 

Here s and C are stretch measure and color measure for the SN 
light curve, Al is a parameter representing some combination of 
SN absolute magnitude M and Hubble constant H 0 . 

For a set of N SNe with correlated errors, the x~ function is 
given by 

x\ N - Am r • C- 1 • Am, (22) 


where Am = ihb - m m0 d is a vector with N components, and nin 
is the rest-frame peak B-band magnitude of the SN. The total 
covariance matrix C can be written as Conley et al. ([2011 1 


C — D s tat + C s tat “ C s 


(23) 


Here D stat denotes the diagonal part of the statistical uncertainty, 
C stat and C sys denote the statistical and systematic covariance 
matrices, respectively. For the details of constructing the covari¬ 
ance matrix C, see Conley et al. ( |2011| >. 

It must be emphasized that, in order to include host-galaxy 
information in the cosmological fits, Conley et al. split the 
SNLS3 sample based on host-galaxy stellar mass at 10 111 /VP,, 


and made M to be different for the two samples Conley et al. 
( |201l| >. So there are two values of M (i.e. All and AH) for the 
SNLS3 data. Moreover, Conley et al. removed M\ and AH 
from cosmology-fits by analytically marginalizing over them 
(for more details, see the appendix C of Conley et al. (|2011[>). 
In the present work, we just follow the recipe of |Con 
( 201 1[ >, and do not treat M as model parameter. 


ey et al. 


2.2.2. CMB Data 


For CMB data, we use the distance priors data extracted from 
Planck first data release (Wang & Wang 2013b) CMB give us 


4 It should be mentioned that, in this paper we just use the purely ge¬ 
ometric measurements of CMB, i.e. the distance prior data. There are 
some other methods of using CMB data. For examples, the observed 
position of the first peak of the CMB anisotropies spectrum can be used 
to perform cosmology-fits (Ca meiro et al.[ 2008, Pig ozzo et al.|20lT| >. 
In addition, the CMB full data can also be used to constrain cosmo¬ 
logical models via the global fit technique. To make a comparison, we 
constrain the parameter spaces of IwCDMl model by using these three 
methods of using CMB data. We find that the differences among the 
fitting results given by different method are very small (For example, 
the differences on DE EoS w are only the order of 1%). In other words. 
Our results are insensitive to the method of using CMB data. This con¬ 
clusion also holds true for other DE models (such as HDE model (jLi et| 
|al.|2013| », showing that the CMB data cannot put strict constraints on 
the properties of DE. Since the main purpose of using CMB data is to 
put strict constraints on fl c0 an d 0 , we think that the use of the Planck 
distance prior is sufficient enough for our work. 


Article number, page 3 of 


12 






















































the comoving distance to the photon-decoupling surface r(z*|^] 
and the comoving sound horizon at photon-decoupling epoch 
r s (z„). Wang and Mukherjee showed that the CMB shift param¬ 
eters ( |Wang & Mukherjee |2007| > 

l a = nr(z,)/r s (z*), 

R = yjci m H 2 r(z*)/c, (24) 

together with cj b = Q.bh 2 , provide an efficient summary of CMB 
data as far as dark energy constraints go. The comoving sound 
horizon is given by ( jWang & Wang|2013b| > 

r" da' 

r s(z) = cH 0 1 I _ (25) 

Jo ^3(1 +R b a') a' 4 E 2 (z') 


where a is the scale factor of the Universe, R b = 
31500Cl b h 2 (T cmh /2J K)- 4 , and T cmh = 2.7255 K. 

Using the Planck+lensing+WP data, the mean values and lcr 
errors of {l a ,R,to b } are obtained (Wang & Wang[2013b|, 


(l a ) = 301.57, <x(/ a ) = 0.18, 

(R) = 1.7407 ,cr(R) = 0.0094, 

(a> b ) = 0.02228, o-(u) b ) = 0.00030. (26) 


Defining p\ = l a (z*), pi = R(z.,), and /?? = cjb, the normalized 


covariance matrix NormCovcMBCp;, Pi) can be written as (Wang 
& Wang|2013bl ) 


1.0000 0.5250 -0.4235 f 

0.5250 1.0000 -0.6925 

, -0.4235 -0.6925 1.0000 , 


(27) 


Then, the covariance matrix for (l a , R, oji,) is given by 


Co v C MB(Pi, Pj) = o-(pt) a-(pj) NormCo v CMB (Pi, Pj), (28) 


where i,j = 1,2,3. The CMB data are included in our analysis 
by adding the following term to they ' 2 function: 


X 2 cmb = A Pi [C° v cmb(P<, Pj)] A Pj . A Pi = Pi - P d i a ' a ’ 


where pf‘ M are the mean values from Eq. (|26 1 . 


(29) 


2.2.3. GC Data 


To improve the cosmological constraints, we also use the GC 
data extracted from SDSS samples. In ( jChuang & Wang |2013| >, 
Chuang and Wang measured the Hubble parameter H(z) and 
the angular-diameter distance D,\ (z) separately |^] by scaling the 
model galaxy two-point correlation function to match the ob¬ 
served galaxy two-point correlation function. Since the scal¬ 
ing is measured by marginalizing over the shape of the model 
correlation function, the measured H(z) and I) A (z) are model- 
independent and can be used to constrain any cosmological 
model ( [Chuang & Wang |2013( l. 

ft should be mentioned that, compared with using H(z) and 
D a (z), using H(z)r s (zj)/c and D A {z)lr s {zd )^\ can give better con¬ 
straints on various cosmological models. In (Hemantha, Wang 


5 In this work, we adopt the result 


oi z » given 




1996). 


^^ngular-diameterdistance^D£(z)^^c^T|C)/(J^^)^here^cjsJhe 


speed of light. 


7 r s (zd) isThe sound horizon at the drag epoch, where r s (z) is given by 

Eq. |25 Zd is given in {Eisenstein & Hu 
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|& Chuang|2014| , using the two-dimensional matter power spec¬ 
trum of SDSS DR7 samples, Hemantha, Wang, and Chuang got 

H(z = 0.35 )r s ( Zd )/c = 0.0431 ± 0.0018, 

D a (z = 0.35)/r s (z d ) = 6.48 + 0.25. (30) 

In a similar work (Wang ||2014) l, using the anisotropic two- 
dimensional galaxy correlation function of SDSS DR9 samples, 
Wang obtained 


H(z = 0.57 )r s (z d )/c = 0.0444 ± 0.0019, 

D a (z - 0.57)/r s (z d ) = 9.01+0.23. (31) 

GC data are included in our analysis by adding Xcc ~ Xcci + 
Tgc 2 ’ z oci = 0.35 and zcci = 0.57, to the x 2 of a given 
model. Note that 


xla = Ac n [CccMi’ qj)] Ac ih Ac n = <n ~ <rf ata > (32) 


where qi = H(z c C i)r s (zd)/c, q 2 = D A (zca) / rAzd), and i - 1,2. 
Based on Refs. (Hemantha, Wang & Chuang [2014] ) and ( jWang | 
|2014j ), we have 


Cgci - 


0.00000324 -0.00010728 

-0.00010728 0.0625 


(33) 


Ccc 2 - 


0.00000361 0.0000176111 

0.0000176111 0.0529 


2.2.4. Totals 2 Function 
Now the totals 2 function is 

X — Xsn + Xcmb + Xgc- 


(34) 


(35) 


We perform an MCMC likelihood analysis using the “Cos- 
moMC” package ( [Lewis & Bridle|2002[ ). 


2.3. Dark Energy Diagnosis and Cosmic Age 


Since previous works pointed out that the differences among the 
cosmological results given by the three different SNLS3 LCF 
are rather small, we need more tools to distinguish the effects of 
different LCF. In the present work, we use the Hubble parameter 
H(z), the deceleration parameter q(z) . the statefinder hierarchy 
5'°(z) and S A \z) (fArabsalmani & Sahni 2009; Zhang, Cui & 
jZhang||2014j ), and the cosmic age r(z) as our diagnosis tools to 
distinguish the effects of different LCF. This subsection consists 
of two parts. Firstly, we introduce the tools of DE diagnosis, 
including the Hubble parameter //(z), the deceleration parameter 


q(z), and the statefinder hierarchy {5 3 ^, Then, we discuss 
the issues about cosmic age, including the high-redshift cosmic 
age test and the fate of the Universe. 

Let us introduce the DE diagnosis tools first. The scale factor 
of the Universe a can be Taylor expanded around today’s cosmic 
age fo as follows: 


CXI 1 

a W = l + Z -rW-'o)]". 

“ n\ 

n= 1 


(36) 


where 


_ a( 0 w 

a(t)H n ’ 


n e N, 


(37) 
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witha(f) (n) = d"a(t)ldt n . The Hubble parameter H(z) contain the 
information of the first derivative of a{t). Based on the Baryon 
Acoustic Oscillations (BAO) measurements from the SDSS data 
release 9 and data release 11, Samushia et al. gave Host = IKz = 
0.57) = 92.4 + 4.5km/s/Mpc Samushia et al. ( 2013| >, while 
Delubacetal. obtained 7 / 2.34 = H(z = 2.34) = 222±7km/s/Mpc 
|Delubac et al.| ( j2014| ). These two H(z) data points will be used 
to compare the theoretical predictions of the wCDM model and 
the three IDE models. In addition, the deceleration parameter q 
is given by 


q = ~A 2 = - 


a 

aH 2 ’ 


(38) 


which contains the information of the second derivatives of a(t). 
For the ACDM model, A 2 |acdm = 1 - §O m , A 3 | A cdm = 1, 

-4 4 1 acdm = 1 ~ The statefinder hierarchy, S is defined as 

|Arabsalmani & Sahni| ( j2009) l: 

S 2 =A 2 + §O m , (39) 

5 3 =A 3 , (40) 

S 4 = A 4 + jfi m , (41) 

The reason for this redefinition is to peg the statefinder at unity 
for ACDM during the cosmic expansion, 


S n Iacdm - 1 ■ 


(42) 


This equation defines a series of null diagnostics for ACDM 
when n > 3. By using this diagnostic, we can easily distin¬ 
guish the ACDM model from other DE models. Because of 
Fl m Iacdm = |(1 + q), when n > 3, statefinder hierarchy can 
be rewritten as: 

Sf = A it (43) 

= A 4 + 3(1 + q), (44) 

where the superscript (1) is to discriminate between S ( J 1 and S „. 
In this paper, we use the statefinder hierarchy S^z) and S ( ^ } (z) 
as our diagnosis tools to distinguish the effects of different LCF. 

Now, let us turn to the issues about cosmic age. The age of 
the Universe at redshift z is given by 


t(z) = 



dz 

(1 +z)H(z)' 


(45) 


In history, the cosmic age problem played an important role in 
the cosmology |Alcaniz & Lima (|1999|l;|Lan et alT] (201 0[>; Ben- 


|galy Jr. et al.| ( |2014] >; |Liu & Zhang| ( |2014) ). Obviously, the Uni¬ 
verse cannot be younger than its constituents. In other words, 
the age of the universe at any high redshift z. cannot be younger 
than its constituents at the same redshift. There are some old 
high redshift objects (OHRO) considered extensively in the lit¬ 
erature. For instance, the 3.5Gyr old galaxy LBDS 53W091 at 
redshift z = 1-55 Dunlop et al. (jl996|, the 4.0Gyr old galaxy 
LBDS 53W069 at redshift z = 1.43 Dunlop (1999), and the 
old quasar APM 08279+5255, whose age is estimated to be 
2.0 Gyr, at redshift z = 3.91 Hasinger, Schartel, & Ko mossa] 
(|2002|). In the literature, the age data of these three OHRO (i.e. 


f 143 = t(z = 1.43) = 4.0Gyr, 1 1.55 = t(z = 1.55) = 3.5 Gyr and 
f 3 .91 = t(z = 3.91) = l.QGyr) have been extensively used to test 


various cosmological models (see e.g. Alcaniz, Lima & Cunha 

(2003 

>; Wei & Zhang (2007); 

Wang & Zhang (2008); Wang, Li 

& Li 

2010); Yan, Liu & Wei 

(2014)). In the present work, we 


will use these three age data to test the wCDM model and the 
three IDE models. 

Another interesting topic is the fate of the Universe. The fu¬ 
ture of the Universe depends on the property of DE. If the Uni¬ 
verse is dominated by a quintessence Caldwell, Dave & Stein- 
|hardt| ( fT998] >; |Zlatev, Wang & SteinhardM 1999) or a cosmologi¬ 
cal constant, the expansion of the Universe will continue forever. 
If the Universe is dominated by a phantom |CaldweH] ( |2002j >; 
|Caldwell, Kamionkowski & Weinberg| ( |2003j >, eventually the re¬ 
pulsive gravity of DE will become large enough to tear apart all 
the structures, and the Universe will finally encounter a dooms¬ 
day, i.e. the so-called “big rip” (BR). Setting x = —ln( 1 + z), we 
can get the time interval between a BR and today 


tBR - to 


f 


dx 

Hix)’ 


(46) 


where tsR denotes the time of BR. It is obvious that, for a Uni¬ 
verse dominated by a quintessence or a cosmological constant, 
this integration is infinity; and for a Universe dominated by a 
phantom, this integration is convergence. We would like to in¬ 
fer, from the current observational data, how far we are from a 
cosmic doomsday in the worst case. So in this work we calculate 
the 2cr lower limits of (br - to for all the four DE models. 


3. Result 

In this section, firstly we show the cosmology-fit results of the 
wCDM model given by various SNLS3 samples without and 
with systematic uncertainties, next we present the fitting results 
of the three IDE models, then we show the cosmic evolutions of 
Hubble parameter H(z), deceleration parameter q(z), statefinder 
hierarchy 5 3* Vz) and .S’^ j (z) according to the fitting results of the 
IDE models, finally we perform the high-redshift cosmic age test 
and discuss the fate of the Universe. Since both three SNLS3 
datasets and a time-varying /3 are considered at the same time, 
all the results of this work are new compared with the previous 
studies. 

3. 1. Cosmology Fits 

In table |T] we present the fitting results of the wCDM model 
given by various SNLS3 samples without and with systematic 
uncertainties. The first row of table [T] shows the fitting results 
for the case of only considering constant f3 and statistical un¬ 
certainties. From this row we find that the best-fit values of 
£I,„o and w given by the Combined sample are in between the 
best-fit results given by the SALT2 sample and by the SiFTO 
sample. This result is consistent with the result of jConley et al.| 
( |2011] >. The second row of table [T] shows the fitting results for 
the case of considering constant /3 and statistical+systematic un¬ 
certainties. From this row we can see that, once the systematic 
uncertainties of SNLS3 samples are taken into account, the best- 
fit values of G„,o and w given by the Combined sample are not 
in between the best-fit results given by the SALT2 sample and 
by the SiFTO sample any more. Therefore, the reason causing 
this strange phenomenon is due to the systematic uncertainties of 
SNLS3 samples. To further study this issue, in the third row of 
table [T] we present the fitting results for the case of considering 
linear [j and statistical+systematic uncertainties. We find that, 
after considering the evolution of / 3 , the best-fit value of w given 
by the Combined sample is in between the results given by the 
SALT2 sample and by the SiFTO sample, while the differences 
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Table 1. Fitting results for the wCDM model given by various SNLS3 samples without and with systematic uncertainties. Both the best-fit values 
and the lcr errors of Q,„o and w are listed. “Stat Only” and “Stat Plus Sys” represent the SN data only including statistical uncertainties and the 
SN data including both statistical and systematic uncertainties, respectively. Moreover, for comparison, both the cases of constant and linear [3 
are taken into account. 
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Table 2. Fitting results for the wCDM model and the three IDE models, where both the best-fit values and the 1 a errors of various pa¬ 
rameters are listed. “Combined”, “SALT2” and “SiFTO” represent the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC and the 
SN(SiFTO)+CMB+GC data, respectively. 
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among the best-fit values of Q m(l given by the three SNLS3 sam¬ 
ples are effectively reduced. These results imply the importance 
of considering /Ts evolution in the cosmology-fits. Therefore, 
from now on we only consider the case of linear / 3. 

In table[2] by adopting a linear /j, we give the fitting results of 
the wCDM model and the three IDE models. From this table we 
see that, for all the four DE models, /? significantly deviates from 
a constant, consistent with the results of | Wang, et al.| ( j2014| ; in 
addition, there is no evidence for the existence of dark sector in¬ 
teraction. Although the best-fit results of w given by the three 
LCF are always less than -1, w = -1 is still consistent with 
the current cosmological observations at 2cr CL. Moreover, we 
check the impacts of different SNLS3 LCF on parameter esti¬ 
mation and find that for all the four DE models: (1) The “Com¬ 
bined” sample always gives the largest w; in addition, the values 
of w given by the “SALT2” and the “SIFTO” sample are close to 
each other. (2) The effects of different LCF on other parameters 
are negligible. It is clear that these results are insensitive to a 
specific DE model. 


In Fig. [T] we plot the probability contours at the lcr and 
2cr CL in the y-w plane, for the three IDE models. The best- 
fit values of {y, w] of the “Combined”, the “SALT2” and the 
“SiFTO” data are marked as a black square, a red triangle and 
a green diamond, respectively. For comparison, the fixed point 
{y, w] = {0,-1} for the ACDM model is also marked as a blue 
round dot. A most obvious feature of this figure is that, for all 
the IDE models, the fixed point {0,-1} of the ACDM model is 
outside the lcr contours given by the three data sets; however, 
the ACDM model is still consistent with the observational data 
at 2cr CL. Moreover, according to the evolution behaviors of 
Pde (see Eqs. 10 12 and p~6]> at z —» -1, we divide these y-w 


planes into two regions: the region above the dividing line de¬ 
notes a quintessence dominated Universe (without big rip), and 
the region below the dividing line represents a phantom domi¬ 
nated Universe (with big rip). We can see that, although all the 
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best-fit points given by the three data sets correspond to a phan¬ 
tom, both phantom, quintessence and cosmological constant are 
consistent with the current cosmological observations at 2cr CL. 
This means that the current observational data are still too lim¬ 
ited to indicate the nature of DE. 


3.2. Hubble Parameter, Deceleration Parameter and 
Statefinder Hierarchy 

The lcr confidence regions of Hubble parameter H(z ) at redshift 
region [0,4] for the wCDM model and the three IDE models 
are plotted in Fig. 2} where the two ll(z) data points, H 0 57 and 
7 / 2 . 34 , are also marked by diamonds with error bars for compar¬ 
ison. We find that the data point Z/ 0.57 can be easily accommo¬ 
dated in the wCDM model and the three IDE models, but the 
data point 7 / 2.34 significantly deviates from the lcr regions of all 
the four DE models. In other words, the measurement of 7 / 2.34 is 
in tension with other cosmological observations and this result is 
consistent with the conclusion of|Sahni, Shafieloo & Starobinsky 
p0l4] l; [Hu et al. ( |2014| l. In addition, the lcr confidence regions 
of I Hz.) given by different LCF are almost overlap; this means 
that using Il(z) diagram is almost impossible to distinguish the 
differences among different SNLS3 LCF. 

In Fig. [3] we plot the lcr confidence regions of deceleration 
parameter q(z) at redshift region [0,4], for the wCDM model and 
the three IDE models. Again, we see that the lcr confidence re¬ 
gions of q(z) given by different LCF are almost overlap. This 
implies that using q(z) diagram also has great difficulty to distin¬ 
guish the differences among different SNLS3 LCF. 

In Fig. J 4 J we plot the lcr confidence regions of S^\z) at 
redshift region [0,4], for the wCDM model and the three IDE 
models. From this figure we see that, the lcr confidence regions 
of Sj'Vz) given by the three SNLS3 LCF are almost overlap at 
high redshift. Although there is an separating trend for .S'^Cz) 
when z —» 0, most parts of these three lcr regions are overlap 
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IwCDMI IwCDM2 




IwCDM3 



V 

Fig. 1. (color online). Probability contours at the lcr and 2cr CL in the y-w plane, for the IwCDMI (upper left panel), the IwCDM2 (upper 
right panel) and the IwCDM3 (lower panel) model. “Combined” (black solid lines), “SALT2” (red dashed lines), and “SiFTO” (green dash-dotted 
lines) denote the results given by the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC, and the SN(SiFTO)+CMB+GC data, respectively. 
Furthermore, the best-fit values of {y, w) of the “Combined”, the “S ALT2” and the “SiFTO” data are marked as a black square, a red triangle and 
a green diamond, respectively. To make a comparison, the fixed point jy, w) = {0, — 1} for the ACDM model are also marked as a blue round dot. 
The gray solid line divides the panel into two regions: the region above the dividing line denotes a quintessence dominated Universe (without big 
rip), and the region below the dividing line represents a phantom dominated Universe (with big rip). 


even at current epoch. This means that, although S^'Vz) is a bet¬ 
ter tool than H(z) and q(z), it still have difficulty to distinguish 
the effects of different SNLS3 LCF. Moreover, our conclusion 
holds true for all the four DE models, showing that this conclu¬ 
sion is insensitive to a specific interaction form. 

In Fig. JtJ we plot the lcr confidence regions of S { ^\z) at 
redshift region [0,4], for the wCDM model and the three IDE 
models. From this figure we see that most of the lcr confidence 
regions of S ^\z) given by the three LCF are overlap. This means 
that the effects of different LCF can not be distinguished by using 


the statefinder S ( ^\z) either. Again, we see that this conclusion 
is insensitive to a specific interaction form. 

In conclusion, we find that the differences given by 
“SALT2”, “SiFTO”, and “Combined” LCF are rather small and 
can not be distinguished by using H(z), q(z), S { ^\z) and S { ^\z). 
This result is quite different from the case of “MLCS2k2” |Jha| 


et 


al. ( |2007j ) and “SALT2” Guy et al. ( 2007| >, where using 
MLCS2k2” and “SALT2” LCF will give completely different 
cosmological constraints for various models Bengochea]( [2011| ); 
Bengochea & De Rossi ( 2014| >. 


3.3. Cosmic Age and Fate of The Universe 

The 2cr confidence regions of cosmic age t(z ) at redshift region 
[0,4] for the wCDM model and the three IDE models are plotted 
in Fig. [ 6 ] where the three t(z) data points, t\. 43 , U .55 and (3 91 , are 
also marked by Squares for comparison. We find that both C .43 


and fi .55 can be easily accommodated in all the four DE models, 
but the position of *3 91 is significantly higher than the 2 cr upper 
bounds of all the four DE models. In other words, the existence 
of the old quasar APM 08279+5255 still can not be explained 
in the frame of IDE model. This result is consistent with the 
conclusions of previous studies Alcaniz, Lima & Cunha (|2003|>; 
Wei & Zhang ( 2007|>; Wang & Zhang| (|2008|l; Wang, Zhang &| 


Xia| ( |2008[ ); |Wang, Li & Li| ( |2010| l; |Yan, Liu & Wei| < ]2014| >. In 
addition, the 2cr regions of t(z ) given by different LCF are almost 
overlap, showing that the impacts of different SNLS3 LCF can 
not be distinguished by using the age data of OHRO. 

We want to infer how far we are from a cosmic doomsday 
in the worst case. So in Fig. [7] we plot the 2 cr lower limits of 
the time interval t - to between a future moment and today, for 
the viCDM model and the three IDE models. Moreover, the 20 - 
lower limit values of tw? - to given by the three LCF are also 
marked on this figure, rl All the evolution curves of t - to tend 
to the corresponding convergence values at -ln( 1 + z) — 20. The 
most important factor of determining Ibr - to is the EoS w. In 
a phantom dominated Universe, a smaller w corresponds to a 
larger increasing rate of p,i e \ this means that all the gravitation¬ 
ally bound structures will be torn apart in a shorter time, and the 
Universe will encounter a cosmic doomsday in a shorter time, 
too. Among the three SNLS3 LCF data, the “Combined” sam¬ 
ple always gives the largest vv, and thus gives the largest value of 


8 The 2cr upper limits of t BR - 1 0 is infinite; in other words, the Universe 
will expand eternally 
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Fig. 2. (color online). The lcr confidence regions of Hubble parameter H{z) at redshift region [0,4], for the wCDM (upper left panel), the 
IwCDMl (upper right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model, where the data points of Ho 57 and 
H z 34 are also marked by diamonds with error bars for comparison. “Combined” (gray filled regions), “SALT2” (blue solid lines), and “SiFTO” 
(purple dashed lines) denote the results given by the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC, and the SN(SiFTO)+CMB+GC 
data, respectively. 


tBR ~ to- In addition, the values of w given by the “SALT2” and 
the “SIFTO” sample are close to each other, and thus the values 
of tBR - to given by these two samples are close to each other. 


4. Summary and Discussion 


As is well known, different LCF will yield different SN sam¬ 
ple. In 2011, based on three different LCF, the SNLS3 group 
Conley et al. (201 ljl released three kinds of SN samples, i.e., 
“SALT2”, “SiFTO” and “Combined”. So far, only the SNLS3 
“Combined” sample is studied extensively, both the “SALT2” 
and the “SiFTO” data sets are seldom taken into account in the 
literature. Moreover, the cosmological consequences given by 
these three SNLS3 LCF are not discussed before. Therefore, the 
impacts of different SNLS3 LCF have not been studied in de¬ 
tail in the past. The main aim of the present work is presenting 
a comprehensive and systematic investigation on the impacts of 
different SNLS3 LCF. 

Since the interaction between different components widely 
exist in nature, and the introduction of a interaction between 
DE and CDM can provide an intriguing mechanism to solve 
the “cosmic coincidence problem” and alleviate the “cosmic age 
problem”, here we adopt the wCDM model with a dark sector 
interaction. To ensure that our study is insensitive to a specific 
dark sector interaction, three kinds of interaction terms are taken 
into account in this work. In addition, to make a comparison, 
we also consider the case of wCDM model without dark sector 
interaction. 

We have used the three SNLS3 data sets, as well as the ob¬ 
servational data from the CMB and the GC, to constrain the pa¬ 
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rameter spaces of the wCDM model and the three IDE models. 
According to the results of cosmology-fits, we have plotted the 
cosmic evolutions of Hubble parameter H(z). deceleration pa¬ 
rameter q(z), statefinder hierarchy S^iz) and S { ^\z), and have 
checked whether or not these DE diagnoses can distinguish the 
differences among the results of different LCF. Furthermore, we 
have performed high-redshift cosmic age test using three OHRO, 
and have explored the fate of the Universe. 

We find that for the wCDM and all the three IDE models: (1) 
the “Combined” sample always gives the largest w; in addition, 
the values of w given by the “SALT2” and the “SIFTO” sample 
are close to each other, (see table [2] and Fig. [TJ; (2) the effects 
of different SNLS3 LCF on other parameters are negligible (see 
table [2). Besides, we find that the ACDM model is inconsistent 
with the three SNLS3 samples at lcr CL, but is still consistent 
with the observational data at 2cr CL. 

Moreover, we find that the impacts of different SNLS3 LCF 
are rather small and can not be distinguished by using H(z ) (see 
Fig. B, q(z )(see Fig. B, Sl^(z), (see Fig. [4J), S^\z) (see FigjiJ, 
and t\z) diagram (seeTig. |6j. This result is quite different from 
the case of “MLCS2k2”|Jhaet al. (2007) and “SALT2”|Guy et al.| 
( |2007| l, where using “MLCS2k2” and “SALT2” LCF will give 
completely different cosmological constraints for various mod¬ 
els Bengochea ( |2011) l; |Bengochea & De Rossi| ( |2014[ l. In addi¬ 
tion, we infer how far we are from a cosmic doomsday in the 
worst case, and find that the “Combined” sample always gives 
the largest 2cr lower limit of - to, while the results given by 
the “SALT2” and the “SiFTO” sample are close to each other 
(see Fig. [7j. 
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Wang el al.: Impacts of different SNLS3 LCF on cosmological consequences of IDE models 





Fig. 3. (color online). The 1 a confidence regions of deceleration parameter q(z) at redshift region [0,4] for the wCDM (upper left 
panel), the IwCDMl (upper right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model. “Combined” (gray 
filled regions), “SALT2” (blue solid lines), and “SiFTO” (purple dashed lines) denote the results given by the SN(Combined)+CMB+GC, the 
SN(SALT2)+CMB+GC, and the SN(SiFTO)+CMB+GC data, respectively. 






Fig. 4. (color online). The lcr confidence regions of S^\z) at redshift region [0,4] for the wCDM (upper left panel), the IwCDMl (upper 
right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model. “Combined” (gray filled regions), “SALT2” (blue 
solid lines), and “SiFTO” (purple dashed lines) denote the results given by the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC, and the 
SN(SiFTO)+CMB+GC data, respectively. 
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Fig. 5. (color online). The lcr confidence regions of S^(z) at redshift region [0,4] for the wCDM (upper left panel), the IwCDMl (upper 
right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model. “Combined” (gray filled regions), “SALT2” (blue 
solid lines), and “SiFTO” (purple dashed lines) denote the results given by the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC, and the 
SN(SiFTO)+CMB+GC data, respectively. 
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Fig. 6. (color online). The 2cr confidence regions of cosmic age t{z) at redshift region [0,4], for the wCDM (upper left panel), the IwCDMl 
(upper right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model. Three t(z) data points, q 43 , t l 55 and t 3 91 , are also 
marked by Squares for comparison. “Combined” (gray filled regions), “SALT2” (blue solid lines), and “SiFTO” (purple dashed lines) denote the 
results given by the SN(Combined)+CMB+GC, the SN(SALT2)+CMB+GC, and the SN(SiFTO)+CMB+GC data, respectively. 
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-ln(1+z) -ln(1+z) 


Fig. 7. (color online). The 2 a lower limits of the time interval t — to between a future moment and today, for the wCDM (upper left 
panel), the IwCDMI (upper right panel), the IwCDM2 (lower left panel) and the IwCDM3 (lower right panel) model. “Combined” (black 
solid lines), “SALT2” (red dashed lines), and “SiFTO” (green dash-dotted lines) denote the results given by the SN(Combined)+CMB+GC, 
the SN(SALT2)+CMB+GC, and the SN(SiFTO)+CMB+GC data, respectively. The 2cr lower bound values of t BR - t 0 given by the three SNLS3 
data are also listed on this figure. 


Since the conclusions listed above hold true for all the three 
IDE models, we can conclude that the impacts of different 
SNLS3 LCF are insensitive to the specific forms of dark sec¬ 
tor interaction. In addition, these conclusions also come into 
existence for the case of the wCDM model. Our method can be 
used to distinguish the differences among various cosmological 
observations (e.g., see Hu et al. (2015 i;jWang et al.](|2016a)i). 


For simplicity, in the present work we only adopt a constant 
w, and do not consider the possible evolution of w. In the liter¬ 
ature, the dynamical evolution of EoS are often explored by as¬ 


suming a specific ansatz for w(z) Chevallier & Polarski (2001); 
Linder P003|l; Gerke, & Efstathiou| (|2002| i; |Wetterich| ( 2004 i; 
Jassal, Bagla, & Padmanabhan (2005 , or by adopting a binned 
parametrization Huterer, & Starkman 2003); Huterer, & Cooray 


(2005 1 ; Huang et al.| ( 2009 1 ; [Wang, Li & Li| ~i 2011 [ i; Li et al. 
(2011 ; |Gong, Gao, & Zhuf ( 2013| l. To further study the impacts 
of various systematic uncertainties of SNe la on parameter esti¬ 
mation, we will extend our investigation to the case of a time- 
varying w in the future. 


In a recent paper |Betoule et al7| ( |2014j ), based on the improved 
SALT2 LCF, Betoule et al presented a latest SN data set (“joint 
light-curve analysis” (JLA) data set), which consists of 740 SNe 
la. Adopting a constant a and a constant /3, Betoule et al. found 
Q,„o = 0.295 + 0.034 for a flat ACDM model; this result is differ¬ 
ent from the result of SNLS3 data, but is consistent with the re¬ 
sults of Planck Ade et al. ( |2014[ >. It would be interesting to apply 
our method to compare the differences between the SNLS3 and 
the JLA sample. These issues will be studied in future works. 
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